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Abstract. We investigate the dependence of galaxy clustering on the galaxy intrinsic luminosity at high redshift, 
using the data from the First Epoch VIMOS-VLT Deep Survey (VVDS). The size (6530 galaxies) and depth 
(7 ab < 24) of the survey allows us to measure the projected two-point correlation function of galaxies, Wp{rp), 
for a set of volume-limited samples up to an effective redshift {z) = 0.9 and median absolute magnitude —19.6 < 

Mb < —21.3. Fitting Wp{rp) with a single power-law model for the real-space correlation function ^(r) = (r/ro) ', 
we measure the relationship of the correlation length ro and the slope 7 with the sample median luminosity for the 
first time at such high redshift. Values from our lower-redshift samples (0.1 < 2 < 0.5) are fully consistent with 
the trend observed by larger local surveys. In our high redshift sample (0.5 < 2 < 1.2), we find that the clustering 
strength suddenly rises around Mb, apparently with a sharper inflection than at low redshifts. Galaxies in the 
faintest sample {{Mb) = —19.6) have a correlation length ro = 2 . 7 jjQ '3 h~^ Mpe, compared to ro = 5.0jl)'g h~^ 

Mpe at {Mb) = —21.3. The slope of the correlation function is observed to correspondingly steepen significantly 
from 7 = l.blo 3 to 7 = 2 . 4 ji;Q’ 2 . This is not observed either by large local surveys or in our lower-redshift samples 
and seems to imply a significant change in the way luminous galaxies trace dark-matter halos at 2 ~ 1 with 
respect to 2 ~ 0. At our effective median redshift 2 ~ 0.9 this corresponds to a strong difference of the relative 
bias, from b/b* < 0.7 for galaxies with L < L* to b/b* ~ 1.4 for galaxies with L > L*. 

Key words, cosmology: deep redshift surveys -large scale structure of the Universe - methods: statistical - galaxies: 
evolution 



1. Introduction 
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ones 
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to be 
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more clustered than faint 
ih, 1988[ JHamilton, 1988 


lovino et ah, 19931 |M 

aurogordato & Lachieze-Rey, 1991 

Benoist et ah, 1996| 

Willmer et ah, 1998| 

Guzzo et 


ah, 2000, . 

IZehavi et ah, 200^ with the difference becoming re¬ 
markable above the characteristic luminosity L* of 
the Schechter luminosity function. This effect is in 
general agreement with predictions from hierarchi- 
cal models of galaxy formation ([White et ah, 1987| 

Valls-Gabaud, Alimi fc Blanchard 198^ Kauffmann et al. 

1997, |BensoireniL7 20011 ), in which bright galaxies are 
expected to occupy more massive dark matter haloes than 
fainter ones and these haloes are more strongly clustered 
than the overall distribution of dark matter ||Kaiser, 1984| 


|Mo fc White, 199'6| |Sheth fc Tormen, 1999 ). If this is 
the case, the difference in clustering between faint and 
bright galaxies should become even more evident at high 
redshifts, where galaxy formation is supposed to be more 
confined to the highest peaks of the density field. 

Understanding the relationship between galaxies and 
dark matter halos is one of the most difficult challenges 
of the theory in predicting the observed clustering 
of galaxies. Over the last few years, halo occupation 
models have provided this connection in a phenomeno¬ 
logical way, allowing one, e.g., to explain the detailed 
shape of the galaxy two-point correlation function 
IZehavi et ah, 20d4| but see also |Guzzo et ah, 19911 ). 
In these models, a statistically motivated recipe to 
describe galaxy formation determines the halo occupation 
distribution (HOD), specifying the probability P(A^|M) 
that a halo of virial mass M contains N galaxies of a 
given type, together with any spatial and velocity biases 
of galaxies | Kauffmann et ah , J^99(^ Benson et ah, 2dMl 


IBerlind et al., 2003r ~ nfeivtso^^^iLi20Mr This term 

(known as the one-halo component of the correlation 
function) governs the behaviour of galaxy correlations on 
small (< 2 h~^ Mpc) scales, while at larger separations 
galaxy correlations are dominated by the gravitational 
clustering of virialized dark matter halos (the two-halo 
component), with essentially no dependence on the more 
complex physics of the sub-dominant baryonic compo¬ 
nent. Given cosmological parameters and a specified 
HOD, therefore, one can calculate any galaxy clustering 
statistic, on any scale (e.g., |Abazajian et ah, 20051 ), either 
by populating the halos from an N-body simulation (e.g., 
|Jing et ah, 1998| |Jlng^t_aI^_2002]) or via analytical pre¬ 
scriptions le.g.. IPeacock fc SmilT^2000l [SeliaU^^000 
[Marinoni et al.,' lOOl IGooray fc Sheth, 2002| 
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|van den Bosch et al., 2003| ). On the other hand, 
as it has been shown USheth fc Tormen, 2004[ 
|Gao, Springel fc White 200^ |Harker et ah, 2005| ), there 
seems to exist a clear relationship between halo formation 
properties and halo clustering properties, which indicates 
that current HOD models may describe galaxy clustering 
at best approximately. Thus, observations of the relative 
clustering of galaxies with different intrinsic luminosities 
provide crucial constraints on HOD models. 

The detailed luminosity dependence of clustering so far 
has been difficult to establish because of the limited dy¬ 
namic range in luminosity for even the largest local galaxy 
redshift surveys (e.g. |Norberg et ah, 2001| ). It is even more 
problematic to study this effect at redshifts significantly 
different to zero. High redshift samples have been too 
small to allow sub-division into luminosity classes. An ad¬ 
ditional complication relates to evolution of the overall 
luminosity function: galaxies become brighter on average 
going back in time, thus comparison of high-redshift mea¬ 
surements to local values requires accurate knowledge of 
the evolution of the global luminosity function. 

The VIMOS-VLT Deep Survey (VVDS) provides us 
with unique information to address these issues in detail. 
A first investigation of how the non-linear bias between 
galaxy and matter evolves with redshift for different lumi¬ 
nosity classes has been presented in Marinoni et al. (2005). 
An analysis of the evolution of clustering of galaxies has 
been presented in Le Fevre et al. (2005a), and the evolu¬ 
tion of the dependence of clustering on spectral types has 
been discussed by Meneux et al. (2005). In this paper we 
use the same VVDS first-epoch data to measure in more 
detail the dependence of galaxy clustering on luminosity 
at (z) ~ 0.9, and compare it to local values from 2dFGRS 
and SDSS. We describe the VVDS catalog and the con¬ 
struction of volume limited samples in Section 2. Section 
3 presents the methods to estimate and retrieve the best- 
fit parameters for the real-space correlation function. We 
present our results on the projected correlation function 
in Section 4, while the comparison to existing local sur¬ 
veys, together with a discussion of the results is given in 
Section 5. 

Throughout this paper we use a Concordance 
Cosmology with 17^ = 0.3 and Da = 0.7. The Hubble con¬ 
stant is normally parameterised via h = iLg/lOO to ease 
comparison with previous works, while a value Hq = 70 
km s“^ Mpc“^ is used when computing absolute magni¬ 
tudes. All correlation length values are quoted in comoving 
coordinates. 

2. The data 

2.1. The VVDS 

The VIMOS-VLT Deep Survey (VVDS, 
|Le Fevre et ah, 2005b| ) studies the evolution of galax¬ 
ies and the large scale structure of the Universe to 
z ~ 2 and higher. The VVDS spectroscopic survey is 
performed with the VIMOS multi-object spectrograph 




























































A. Polio et al.: The VIMOS-VLT Deep Survey 


3 


-24 


-22 

<-20 

U 

w 

> 

CQ 

-18 


-16 


-14 


_1_^_1_ 

1 i 1 

1 ^ ^ \ 1 ^ 

__J_^_1_ 

- 


• ,f. '■ - ■■•' ] f. 


_ 


S.b';:.' iGrL;. ■i.-. 





- 

- 



- 



- 






-■ 



1 1 

-■ 


, , , , 1 

1 1 


0 0.5 1 1.5 


z 


Fig. 1. Distribution of magnitudes of VVDS galaxies with 
redshift. Solid vertical lines show the boundaries of our 
two redshift slices. Solid horizontal lines show from which 
magnitudes {Mb ~ 19.5) our sub-samples are volume lim¬ 
ited. 


at the European Southern Observatory Very Large 
Telescope ( |Le Fevre et ah, 2003| l and complemented with 
multi-color BVRI imaging data obtained at the CFHT 
telescope ( [McCracken et ah, 2003||Le Fevre et ah, 2004| ). 
The first-epoch VVDS data consist of more than 11000 
spectra obtained in two VVDS-Deep fields, covering 0.61 
square degrees. 

For this analysis, we use the data from the F02 “Deep” 
field, which is a purely magnitude limited survey to Iab = 
24 currently covering an area of 0.49 square degrees. 
Considering only galaxies with secure (> 80% confidence) 
redshifts, this sample includes 6530 galaxies. Details about 
observations, data reduction, redshift measurement and 
quality assessment can be found in Le Fevre et al. (2005b). 

2.2. Luminosity-limited sub-samples 

To measure the dependence of clustering on galaxy lu¬ 
minosity, we select two redshift slices to isolate compa¬ 
rable intervals of cosmic time. By choosing the intervals 
z S [0.1,0.5] and z G [0.5,1.2], we obtain two samples 
covering approximately 3.5 Gy of proper time, in the 
adopted cosmology. The high-redshift slice is compared 
to the low-redshift sample and to local samples like the 
2dF Galaxy Redshift Survey (2dFGRS) and the Sloan 
Digital Sky Survey (SDSS). The low-redshift sample is also 
directly compared to existing local estimates, although 
the small volume and larger redshift limit do not allow 
a one-to-one match. Within each slice, we build a series 


of volume-limited (where statistically possible) or quasi- 
volume-limited sub-samples chosen to contain a statisti¬ 
cally reasonable number of galaxies. 

Absolute magnitudes for VVDS galaxies have been es¬ 
timated by computing the appropriate K-correction via 
a SED fitting technique to the observed multi-band pho¬ 
tometry (see jllbert et ah, 2005| and [Franzetti et al., 2005| 
for details) and to ease comparison to previous work are 
reported here in the VEGA system. 

Due to the apparent magnitude limit of the survey, as 
shown in Figure ^ within the high-redshift slice a true 
volume-limited sample can be defined only for Mb 

— 19.5. Conversely, in the low-redshift slice we can define 
volume-limited sub-samples for Mb ^ —17. Because of 
the smaller volume probed we do not have a large num¬ 
ber of bright galaxies in this latter sample. Therefore, 
as seen from Figure ^ we cannot expect statistically 
sound measurements for galaxies closer than z ~ 0.5 
with Mb —19.5. All details on the selected VVDS sub¬ 
samples are listed in Tables 1 and 2. For each sub-sample, 
the columns give its absolute magnitude limit in the B 
band (1); its median Mb magnitude (2); the Schechter 
characteristic magnitude in that redshift range (3); the dif¬ 
ference between the median and characteristic magnitudes 
(4); number of galaxies (5); median redshift (6); measured 
correlation length and the slope of the correlation function 
(7 and 8). 

In the following sections we use - at different redshifts 

- the characteristic magnitude of the Schechter luminos¬ 
ity function in the B band, Mg, as a reference value. 
Normalizing our median absolute magnitude values at 
each redshift to the corresponding value of Mg provides 
a way to take into account the mean brightening of galax¬ 
ies due to evolution, when comparing samples at different 
epochs. The values ot Mg at each redshift are those esti¬ 
mated from these same data by Ilbert et al. (2005) in the 
AB system, converted here into the VEGA scale. 


3. Estimating the real-space correlation function 

We summarize here the methods applied to derive the real- 
space correlation function and its parameters, described 
extensively in Polio et al. (2005). The galaxy real-space 
correlation length rg and slope 7 from the VVDS-F02 data 
are measured from the projection of the redshift-space cor¬ 
relation function ^(rp,7r), estimated using the Tandy & 
Szalay estimator. 


_ Nr{Nr-1) GG(rp,7r) 
— NciNc-l) RRirj,,Tv) 
Nr- 1 GR{rp,TT} , 
Ng RR{rp,Tr) W 


( 1 ) 


where Nq and Nb are the total numbers of objects in 
the galaxy sample and in a catalog of random points dis¬ 
tributed within the same survey volume and with the 
same redshift distribution and angular selection biases; 
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Table 1. Properties of the VVDS luminosity sub-samples in the range 0.1 < z < 0.5. 


Mb range 

j^median 

M*b 


Ngal 

^median 

ro 

7 

< -16 

-17.99 

-19.97 

1.98 

1330 

0.332 

9 rt:: + U.b 4 

^•^^-0.54 

“T’^gTirT^ 

i.uy_Q 

Volume limited: 

< -17.0 

-18.30 

-19.97 

1.67 

1089 

0.363 

2 07 +^-^^ 

' -0.46 

1 '79+U.ib 

-L. 1 ^_o.l2 

< -17.5 

-18.53 

-19.97 

1.44 

883 

0.373 

Q 1 '7+0.62 

0.1 { _o .62 

-1-00-0.11 

< -18.0 

-18.98 

-19.97 

0.99 
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0.381 

q 9r+0.90 
O.zo_o.92 

1 '79+0.22 

1 - <+-0.18 

< -18.5 

-19.34 

-19.97 

0.63 

475 

0.388 

3.47t°;®° 

i.OO_Q 20 

< -19.0 

-19.70 

-19.97 

0.27 

318 

0.381 

4.25tlf4 

-*-• '*^-0.26 

< -19.5 

-20.00 

-19.97 

-0.03 

201 

0.391 

Q 

0.00 5 26 

1.0U Q 4Q 

Table 2. As Table 1, but for 0.5 < z 

< 1.2 





Mb range 

ji^median 

M*b 
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^median 

ro 
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< -17 

-19.65 
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1.11 

4283 

0.808 

‘ ^-0.27 

1 5 Q+u.uy 

i-Oy_Q 07 
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3736 

0.856 

z,.oy_0.27 

1 p;4+0-08 
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3272 
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z,.yu_Q 35 
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i-oz,o 08 

Volume limited: 
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0.53 

2407 

0.899 

9 Qq + U.iii5 
Z,.yo_Q 35 

1 CQ+U.i:^! 

i-oy_o 09 
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-20.58 

-20.76 

0.18 

1530 

0.914 

0 ^'t+0.46 
'^•^'-0.43 

1.84l°;l^ 

< -20.5 

-20.92 

-20.76 

-0.16 

865 

0.913 

4 77+0-61 
'-0.61 

2 . 00 t°; 1^2 

< -21 

-21.30 

-20.76 

-0.54 

368 

0.920 

O.Ui-i 65 

9 qQ+0.36 
+ -O 0 o. 2 i 


GG{rp,TT) is the number of independent galaxy-galaxy 
pairs with separation perpendicular to the line of sight 
between tt and tt + dn and separation parallel to the line 
of sight between Vp and Vp + drp\ RR(rp, tt) is the number 
of independent random-random pairs within the same in¬ 
terval of separations and GR{rp, tt) represents the number 
of galaxy-random pairs. A total of ~ 40000 random points 
has been used in each computation. 

To derive the real-space correlation function and avoid 
the effect of peculiar velocities which distort the redshift 
space statistics, we integrate ^(r^, tt) along the line of sight 
(see |Davis &: Peebles, 1983| ), computing the projected cor¬ 
relation function. 


pOC 

o{rp) = 2 / ^(rp,7r)d7i 

Jo 


= 2 e 




dy, 


( 2 ) 


where ^ is the real space two-point correlation function 
evaluated at the specific separation r = Jvp + y'^. In prac¬ 


tice the upper integration limit has to be chosen finite, to 
avoid adding noise to the result. After a set of experiments 
we chose its optimal value as 20 h~^ Mpc. If ^(r) is well 
described by a power law, ^(r) = (r/ro)“^, the integral 
can be computed analytically, giving 


Wpirp) =rp[ — 


ro 


rG)r(V) 

r(i) 


(3) 


where P is Euler’s Gamma function. Fitting Wp{rp) for 
separations < 10 h~^ Mpc using the procedure described 
in detail in Polio et al. (2005) provides a best-fitting value 
of ro and 7 for each volume-limited sub-sample. 

The estimate of errors has been performed primar¬ 
ily using a bootstrap resampling of the data. However, 


our detailed error analysis in Polio et al. (2005), based 
on 50 VVDS mock surveys from the GallCS simulations 
HBlaizot et al., 20051 ), indicates that bootstrap errors tend 
to be an underestimate of the true ensemble errors. The 
difficulty we encounter is that luminosity-selected sub¬ 
samples of our mock surveys do not show the same prop¬ 
erties, in terms both of total numbers and scaling of the 
intrinsic clustering with luminosity, as the real VVDS. We 
have no guarantee, therefore, that for this specific appli¬ 
cation the variance among the mock samples represent a 
realistic estimate of the errors in the real data. For this 
reason, we have decided here to estimate error bars on 
Wp{rp) using the bootstrap technique, but correct these 
empirically to include the contribution of cosmic vari¬ 
ance. An average value for this correction has been esti¬ 
mated from the direct comparison of the errors computed 
in both ways for 50 whole mock samples ) |Meneux, 2005[ 
[Polio et al., 200^ . The overall mean effect is to increase 
the size of the error bars on Wp by a factor of ~ 2 for 
the low-redshift samples, and by a factor of ~ 1.3 for the 
high-redshift samples. We applied this correction to all 
our bootstrap estimates of Wp{rp). We have checked that 
our conclusions are robust to the details of this correction: 
even doubling the error bars in both redshift ranges, the 
trends in the values of correlation length and slope that 
we find at z ^ 1 do not change and remain significant. 


4. Results 

Figures and 13 show the results of the power-law fits of 
the projected correlation functions Wp{rp) and the corre¬ 
sponding ro and 7 error contours, for the low- and high- 
redshift samples. The fitting has been performed taking 






























A. Polio et al.: The VIMOS-VLT Deep Survey 


5 



0.1 1 10 2 3 4 5 6 7 


2.5 
2 

1.5 
1 

2.5 
2 

1.5 
1 

2.5 
2 

1.5 
1 

2.5 
2 

1.5 
1 

2.5 
2 

1.5 
1 

2.5 
2 

1.5 
1 

2.5 
2 

1.5 
1 


hrp (Mpc) 


hr„ (Mpc) 


c^ 














Fig. 2. Projected correlation functions Wp{rp) and the associated tq and 7 error countours for the seven volume limited 
catalogs at z < 0.5. 


into account the full covariance matrix of the data, as de¬ 
scribed in Polio et al. (2005). 

4.1. The correlation length tq 

A comparison of Figures [21 and |21 shows qualitatively that 
for the high-redshift samples both the correlation length 
and the slope of the correlation function change with the 
sample luminosity. Note that we do not perform here any 
analysis of the detailed shape of Wp(rp), its evolution and 
its implications for halo occupation models, but limit our¬ 
selves to the simple and robust fit of Wp{rp) with a single 
power law. 

In the left panel of Figure 0 we compare our 
measurements of tq from the low-redshift VVDS sam¬ 
ples to the similar measurements from the 2dFGRS 
dNorberg et ah, 2001| | and the SDSS ( [Zehavi et al., 2005| |. 
A one-to-one comparison is not possible, since the large 
size of local surveys allowed for the measurements in a 
series of disjoint volume-limited surveys with small mag¬ 


nitude intervals, i.e. < L < L2, with LI and L2 being 
their limiting luminosities. In our case we are forced to use 
integral measurements, i.e. samples with L < L2. Still, 
different sub-samples are dominated by galaxies with a 
specific characteristic luminosity that we characterize by 
computing the median absolute magnitude within each 
sample (reported in Tables 1 and 2). In addition, given the 
smaller volume of our low-redshift samples, it is clear that 
rare luminous objects will be under-represented. Finally, 
our low-redshift samples extend to z = 0.5 and thus may 
be regarded as a subsequent redshift bin after 2dFGRS 
and SDSS. Gonsidering these intrinsic limitations , the 
overall trend is consistent with that observed by 2dFGRS 
and SDSS, although our measurements are systematically 
lower. 

With this consistency in mind we analyze our VVDS 
high-redshift (z ~ 0.9) samples and compare it to the lo¬ 
cal values from 2dFGRS and SDSS. This is done in the 
right panel of Figure 01 Galaxies fainter than Mg at high 
redshift are significantly less clustered than their counter- 
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Fig. 3. Projected correlation functions Wp(rp) and the associated rg and 7 error countours for the seven high-redshift 
(0.5 < z < 1.2) samples. 


parts in the present-day Universe, with rg = 2.75 ± 0.27 
h~^ Mpc. At the same time, the clustering strength of 
galaxies brighter than is comparable to that observed 
locally with a correlation length up to rg = 4.77 ± 0.61 
h~^ Mpc. We therefore observe that at redshift z ~ 0.9, as 
luminosity increases above L*, the clustering length sud¬ 
denly rises to values comparable to those observed locally 
for galaxies with similar Mb — M*. 


4.2. The correlation function slope 7 

As we show in Figure 0 the dependence of the correla¬ 
tion function slope 7 on the galaxy intrinsic luminosity 
in the high-redshift sample differs strongly from the lo¬ 
cal measurements. Locally, 7 has a remarkably constant 
value, with 7 ~ 1.7 measured both by the much larger 
local surveys and - in a very consistent way - by VVDS 
in our closer redhift bin z ~ 0.4. Conversely, at high z, 
f (r) clearly steepens with increasing luminosity for all sub¬ 
samples brighter than Mb — Mg -\- 0.5. We find that for 


Mb — Mb > the slope of the correlation function is 
consistent with 7 = 1.55 ± 0.08, while for Mb — Mg < 0 
the slope goes up to 7 = 2.38ta2i- 


4.3. The relative bias of different luminosity classes 

To interpret our results and to compare them to local sur¬ 
veys, we compute the relative bias parameter, b/b*, which 
gives the amplitude of the correlation function relative to 
that of L* galaxies. Consistently with the 2dFGRS anal¬ 
ysis (e.g. [Norberg et ah, 2002| ) we define the relative bias 
of the generic i-th sample with a given median luminosity 
L, with respect to that corresponding to L*, as 


b* 





(4) 


and estimate it at the fixed r = 1 h~^ Mpc scale (see also 
Meneux et al. 2006 for a slightly different definition). 

To apply this formula to the low-redshift samples, we 
need to estimate the values of rg and 7 for Mg galaxies. 
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values from 2dFGRS and SDSS compared to the VVDS measurements at (z) ~ 0.9. 
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Fig. 5. Measured values of the slope of ^(r), 7 , as a function of the median luminosity of galaxies at low redshift 
(left panel) and at z ~ 0.9 (right panel). While locally and up to z ~ 0.5 7 remains practically constant (with the 
best rms fit 7 = 1.73 — 0.02{Mb — M’^), marked as a solid line), in the distant Universe it shows a clear scaling 
with luminosity for galaxies brighter than + 0.5: the spatial correlation function becomes steeper for galaxies of 
increasing luminosity. 
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that we obtain by a linear fit to the observed VVDS data 
in the left panel of Figure 01 (avoiding the very uncertain 
value measured at M^). We then plot the corresponding 
values of b/b* for our samples in Figure El In this figure 
we also plot the 2dFGRS and SDSS data, together with 


the analytic fitting relations provided for them, b/b* = 
0.85 + 0.16L/L* for the 2dFGRS | |Norberg et al., 2002| ) 
and b/b* = 0.85 + 0.15 L/L*-0.04(M-M*) for the SDSS 
IITegmark et al., 2004|). Given the error bars, we can say 
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Fig. 6. The variation of the relative bias in the VVDS sub-samples (filled circles) as a function of luminosity using the 
clustering of L* galaxies as a reference point, compared to the 2dFGRS (open triangles) and SDSS (open squares) local 
results. The dashed and dotted lines show, respectively, the best fits made for the SDDS and 2dFGRS measurements. 


that the low-redshift VVDS results are in good agreement 
with both the SDSS and the 2dFGRS fits. 

In the high-redshift samples the situation is clearly dif¬ 
ferent. As we can see from the right panel of Fig. |H1 the 
value of b/b* rises steeply from low values b/b* ~ 0.5 at 
low luminosities to values that, statistically, are consis¬ 
tent with those of the local samples for L > L* galaxies, 
b/b* ~ 1 — 1.4. At the same time the difference on rel¬ 
ative bias between of galaxies fainter and brighter than 
L* becomes very large. This appears to be an indication 
that going back in time the bias contrast of luminous 
galaxies to the rest of the population becomes stronger 
and is consistent with the fact that fainter galaxies are 
found to be significantly less biased tracers of the mass 
than the L* galaxies even at this relatively high redshift 
) |Marinoni et al., 2005| |. 

5. Summary and discussion 

The projected correlation functions that we have mea¬ 
sured from our set of volume limited sub-samples of the 
VVDS are in general fairly well fitted by a single power- 
law in the range 0.1 < r/h~^Mpc < 10, both for the 
low-redshift and high-redshift samples. This allows us to 
use variations in ro and 7 to characterize the global de¬ 
pendence of clustering on luminosity at high redshift and 
compare it to similar low-redshift results. The observed 
behaviour has strong implications for HOD models, as it 
directly impacts on any recipe for populating dark mat¬ 
ter halos at high redshift. Deviations from the power- 
law shape, although extremely interesting for constraining 
HOD models (e.g. |Zehavi et al., 2005| | are not analyzed in 
this paper and will be the subject of future work. 


We observe that for median redshifts z ~ 0.9 (0.5 < 
z < 1 . 2 ) the clustering length has a low, nearly constant 
value ro — 2.9 h~^ Mpc for luminosities Mb < M^, rising 
suddenly for Mb > M/^ and approaching values ro — 5 
h~^ Mpc similar to those of local galaxies with compara¬ 
ble luminosity (relatively to the characteristic value M^, 
[Norberg et al., 2001| ). This behaviour is consistent with 
general predictions by hierarchical models of galaxy for¬ 
mation I jBenson et ah, 2001| ) , where luminous galaxies are 
more confined to the peaks of the large-scale density field, 
going back in redshift, simply due to the higher bias of 
the parent halos. 

Another important result of this work is the clear de¬ 
tection of a systematic steepening of the high-redshift cor¬ 
relation function for absolute magnitudes brighter than 
^M*s + 0.5 (Figure EJ. This kind of behaviour is in gen¬ 
eral not seen either in our closer z ~ 0.4 sample or in 
the large local surveys, although Zehavi et al. (2005) do 
detect an increase of 7 in their most luminous volume- 
selected sample. A similar trend with increasing (UV) lu¬ 
minosity has been recently observed for a population of 
Ly-break galaxies at z > 4 in the Subaru Deep Field 
I jKashikawa et al., 2005| ). These authors are able to repro¬ 
duce their observed relationship with a HOD model in 
which they introduce multiple LBGs into massive dark 
matter haloes. This amplifies the clustering strength at 
small scales, steepening the correlation function. A simi¬ 
lar interpretation could be applied to our data. Finally, the 
observed relative bias of galaxies at high redshift provides 
evidence for a clear difference in the clustering properties 
of galaxies fainter or brighter than the characteristic lu¬ 
minosity: sub-samples with Mb ^ Mg behave in a way 
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that is very similar to local samples while the relative bias 
of samples with Mb ^ remains significantly lower. 

Results presented in this paper show that there is 
a significant redshift evolution of the luminosity depen¬ 
dence of both the normalization and slope parameter of 
the galaxy correlation function. This specific observation 
can provide an important test of galaxy formation mod¬ 
els, constraining in particular the multiplicity of luminous 
galaxies within massive halos at z = 1. 

NOTE ADDE D IN PROOF: In a parallel paper 
Coil et al. (2006)] perform a similar measurement at z ~ 1 
using the DEEP-2 survey. Although they explore a nar¬ 
rower range in median luminosities, they also detect a first 
hint of the steepening of Wp{rp) above and a rise of 
the correlation length with luminosity. Considering cosmic 
variance and the different selection function (unlike the 
VVDS, DEEP-2 is not a purely magnitude-limited sur¬ 
vey) , the overall results from these two data sets are thus 
in good agreement. 
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